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Abstract 

We investigate the B s — ► fi + fi~ decay in the presence of a light stabilized radion in Randall- 
Sundrum model. The branching ratio BR{B S — * [i + fi~) in the standard model is found to be 
3.17 x 10~ 9 (two order smaller than the experimental upper bound) and raises the question whether 
some new physics can play a crucial role or not. We found that for a reasonable range of parameters 
(i.e. radion mass and radion vev (</>)), the above deficit between the experimental bound and 
the standard model result can well be accomodated. Using this upper bound on BR(B S — * 
we obtain the lower bound on ((f)). 
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1 Introduction 



The standard model(SM) despite it's enormous experimental success, succumbs some serious draw- 
backs as far as our understandings of the flavour structure of the quarks and leptons are concerned. 
The CKM matrix, telling a lot about the mixing and CP-violation in the quark sector, lacks any 
dynamical mechanism in it's origin. Simlilarly, the neutrino masses and mixing which is absent in 
the SM, requires the existence of beyond SM physics. An urge for going beyond the SM (by invoking 
new physics(NP)) has become a driving force behind most of the present phenomenological studies. 
Ideas like supersymmetry (with or without i?-parity), technicolor, extra-dimensional field theory, 
split fermion as a candidate of NP draw a lot of attention among the physics community and the 
result is that the LANL archive is flourished by thousands of paper in each year. Among several 
ingenious ideas the notion of extra spatial dimension(s) ([1], [2]) have become very popular as they 
partially resolve hierarchy puzzle. Among these the Randall-Sundrum(RS) model (of warped extra 
spatial dimension ) is particularly interesting from the phenomenological point of view. This model 
views the world as 5-dimensional and it's fifth spatial dimension is S* 1 jZi orbifold. The metric 
describing such a world can be written as 

ds 2 = tfrj^dxV - R 2 c d6 2 . (1) 

The pre-factor VL 2 = e ~ 2kR c\S\ appearing in above is called the warp factor. In £l 2 , k stands for 
the bulk curvature constant and R c for the size of the extra dimension. The angular variable 9 
parametrizes the fifth dimension (xs = R c 9). The model is constructed out of two D% branes 
located at the two orbifold fixed points 9 = and 9 = it, respectively. The brane located at 9 = 
(where gravity peaks) is known as the Planck brane, whereas the brane located at 9 = it (on which 
the SM fields resides and the gravity becomes weak) is called the TeV brane. 

The radius R c can be related to the vacuum expectation value (vev) of some modulus field T(x) 
which corresponds to the fluctuations of the metric over the background geometry (given by R c ). 
Replacing R c by T(x) we can rewrite the RS metric at the orbifold point 9 = -it as 

ds 2 = g™dxV, (2) 

where g™J = e~ 27rkT ^r]^ v = (^y^j V^u- Here f 2 = 24 ^5 5 , where M5 is the 5-dimensional Planck 
scale [2] and 4>{x) = f e ~' KkT { x ) m The scalar field 4>(x) (i.e. 4>(x) = 4>(x) — (</>)) is known as the 
radion field. In the minimal version of the RS model there is no potential which can stabilize the 
modulus field T(x) (and thus the radion (f>(x)). However in a pioneering work, Goldberger and 
Wise [3] were able to generate a potential of this modulus field (by adding an extra massive bulk 
scalar) field which has the correct minima satisfying kR c ~ 11 — 12, a necessary condition for the 
hierarchy resolution. 

In this non-minimal RS model (RS model together with the Goldberger and Wise mechanism), 
the stabilzed radion can be lighter than the other low-lying gravitonic degrees of freedom and will 
reveal itself first either in the direct collider search or indirectly through the precission measurement. 
Studies based on the observable consequences of radion are available in the literature ([4]). 

However, the impact of radion in heavy B meson decay, particularly in the rare decay mode, 
is not fully explored and the present work is an effort in that direction. We have investigated the 
rare B s — > mode in the light of a stabilized radion. The branching ratio BR(B S — > [i + 
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within the SM is found to be two order smaller than the experimental upper bound 1.0 x 10~ 7 [5]. 
It is worthwhile to see whether this deficit can be accommodated within this brane world model. 

The organization of the paper is as follows. In Section 2 we obtain the BR(B S — > 
within the standard model and in the presence of a light stabilized radion. We discuss several input 
parameters: CKM matrix elements, quark and lepton masses, effective wilson coefficients, decay 
constants in Section 3. Section 4 is fully devoted to the numerical analysis. Here we show how the 
upper bound of BR(B S — ► fj, + fi~) can give some lower bound on radion vev (<p) for a relatively 
light radion. We summarize and conclude in Section 5. 



2 The exclusive rare B s —> decay 
2.1 Standard Model(SM) 

The B s (pb) — > /U + (piV~ fe) (partonically described as the bs — > FCNC transition) decay 

within the SM is found to be two order smaller than the experimental upper bound. This channel 
can be useful to test certain class of New Physics, particularly the notion of warped geometry by 
probing it's moduli (radion). 

To begin with let us recall the result of the B s — ► /U + /i~ decay in the standard model [6]. The 
QCD improved effective Hamiltonian Ti. e g describing such a AB = 1 transition can be written as 

H cS = ^ [v tb V* (c e 7 ff 7 + c e 9 ff 9 + doOio) ] , (3) 

where the operators Oi{ i = 9,10) (semileptonic operators involving electro- weak (j,Z) penguin 
and box diagram) and O7 (magnetic penguin) [6, 7] are given by, 

o 9 = -(s^p L b)(ji- 7 y+), 

7T 



(X 

Ova = -tWi* p i>V)(P~'f'%li > 

IT 



OL 

7 = -{sa^PRb) 

7T 



-21171), 

q 2 



(/z-tV)- (4) 



Here a^ v = | [7^, 7^], the chiral-projection operators Pr^l = |(1 i 7s)- ot = f - is the QED fine 
structure constant. In above q is the momentum transferred to the lepton pair, m&, the 6-quark 
mass. Since we are considering the B s meson, the matrix element of H e s is to be taken between 
the vacuum and \B S ) state. Defining as the decay constant of the B s meson, we find [7] 

(0\sn 5 b\B s (p B ))=if^ B , (5) 

— % f — 7TI— 

(Ols^blBsipB)) = B ° B > , (6) 
m b + m s 

and 

(0\sa^P R b\B s ( PB )) = 0. (7) 
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Since q = p\ + p 2 , the Cg term in Eqn.(3) gives zero contribution on contraction with the lepton 
bilinear, of^ gives zero by Eqn.(7) and the remaining c\q gives a contribution given by 

iM SM (B s ( PB ) n + {pi)p~{p 2 )) = iM l0 (B s ( PB ) fi + (pi)p-(p 2 )) 

= i^V tb V t * s c 10 (^( Pl )^-(p 2 )\O l0 \B s (p B )}, (8) 

where 0\q is given in Eq. (4). Accordingly the decay amplitude square is given by 

|M 5M | 2 = Prop- (9) 
For the explicit expression of |Mio| 2 see Appendix A. 



2.2 New Physics(NP) (moduli/radion contribution) 

The light stabilized radion <j){x) can potentially be significant in the B s — > p + p~ decay. Since 
gravitational interaction conserves flavour, radion can not cause flavour-changing neutral current 
(FCNC) transition b — > s at the tree level and it can happen only at the loop level which we will see 
shortly. The radion interaction with the SM fields (which lives on the TeV brane) is governed by 
the 4 dimensional general coordinate invariance. It couples to the trace of the energy-momentum 
tensor T^(SM) of the SM fields and is given by 

C mt = j^TH(SM). (10) 



Here ((ft) is the radion vev and Tff(SM) is 



2M^W+W-» - m 2 z Z^ 



+ (2m 2 h h 2 -dphdPh) + ... (11) 



where D^ip = (d^ — ig^-A^ — ig'^-B^)ip (r a 's being the Pauli matrices, g, g' being SU(2)l and 
U(1)y gauge coupling constants and Yyy, the weak hypercharge). Clearly Eq.(ll) manifests that 
radion intercation with matter conserves flavour at the tree level. The Feynman rules comprising 
radion interactions with the SM fields are listed in Appendix B.l (see Figure 1). For vertices involv- 
ing radion-goldstone-goldstone, radion-goldstone-W boson, radion-fermion-antifermion-goldstone, 
radion-fermion-antifermion-W boson, we read the Feynman rule for vertices [8] without the radion 
and then just multiply the factor 1/(0) to get the desired Feynman rule. For other vertices, the 
Feynman rules are directly read from Eq. (10). Now a b quark can decay to a s quark together 
with a ultra-light radion(0, dropping the hat from cj)(x) from now and onwards) followed by the 
<p — > / + /~ (I = e, fi) decay. Since radion coupling to a pair of fermion (on-shell) is proportional to 
the fermion mass and m e = x 10~ 3 , we do not consider the (f> — > e + e~ decay channel (and thus 
B s — > e + e~ decay mode) in our analysis, as it will not lead to useful bound on (cp). At this point we 
should note that a virtual radion of heavier mass (say about few hundred GeV), might appears to 
be an interesting option. However, the possibility of FCNC transition (like b — > s transition) with 
the radion, the real gravi-scalar, is ruled out. Note that an ultra-light radion of mass about 1 ~ 4 
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GeV is not ruled out as far as the neutrino oscillaton inside supernova is concerned (see Ref.[9]). 
In Ref. [9], the authors showed that the radion exchange (potential) between the neutrino and 
the supernova matter will not effect the neutrino oscillation inside the supernova, yielding a lower 
bound on about 1 GeV. Another thing which requires attention is the possibility of having the 
radion (j)(x) mixing with the higgs h(x) field. A term like 



S 



with fl+ = (0 (v + h{x))/y/2^ ( with v = 247 GeV, the EW vev), can cause radion-higgs mixing 
[10]. However, here we confine ourselves in no-mixing scenario ( i.e. we set £, the radion-higgs 
mixing parameter equal to zero). Analysis in a general radion-higgs mixing scenario is underway. 
The Feynman diagrams contributing to the b — > scp FCNC transition are shown in Figure 2 (see 
Appendix B.2). 

It is now straightforward to evaluate those diagrams and the effective operator parametrizing 
the radion contribution to b — > s<f> process can be written as 



On = Vckm -/= £0 [m b sP R b + m s sP L b], 



(12) 



where Vckm = V^V ib (i = u, c, t), G F {~- 
factor takes the form 



8M^ )' ^ ne ^ erm i c constant. C<j>, the effective loop integral 



2tt 2 ((/») 



dx 



4 Log 



TTL; 



4 + x-^- + 
A 2 



(1 



x)M? 



w 



A 2 



(13) 



in the t'Hooft-Feynman gauge ^ = 1 and A 2 = xmf + (1 — x)Myy. In Eq. (13) i runs over u, c and t 
quarks, but in our analysis we set % = t (since mt 3> m u , rrid). In evaluating we assume that there 
is no external momentum flow (actually the external momenta is much smaller than the internal 
masses Mw and mt and hence setting the external momentum to zero is a good approximation). 
Note that, the Figures 2(b,e,i,j) contribute in b — > scf) vertex (Eq. (13)), while the rest are not. We 
use the cut-off regularizing technique in regularizing the loop integral with the UV cut-off being 
A = 4n{4>) (follows from naive-dimensional analysis). 




-> 




Figure 3. Radion contribution to B s (ps 
the semileptonic process b(pb)s(p s ) 
corresponds to the effective b — s — 



• n + (pi)n {p2)decay. Momentum conservation (for 
■ fi + (pi) ^ (p 2 ) ) reads Pb = Pb + Ps = Pi +P2- The red blob 
vertex. 
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In finding the NP (radion) contribution to B s — > p> + p>~ decay, we neglect the m s dependent 
term in Eq.(12) (since m s <C nib). It is now straightforward to include the NP contribution to the 
B s — > p + p~ decay amplitude. The relevant new physics Feynman digram is shown in Figure 3 and 
the squared-amplitude \M$m + -^fiVp| 2 (-Bs — ► H + is given by 



|M 5M + M NP \ 2 = |M 10 | 2 + |M WP | 2 + 2Re(Ml M NP ). (14) 

The direct and interference terms of the amplitude (Eq.(14)) is presented in the Appendix A and 
B.3. Note that, in Figure 3 the radion 0, produced in the annihilation of b and s quarks, decays 
into p^ pT pairs. We are working in the Breit-Wigner approximation and in this approximation 
the radion propagator takes the form 

1 



(15) 



(q 2 - rrify + im^T^ 

Here q 2 = (pb +p s ) 2 = (pi +P2) 2 and is the total width of the radion arising from it's decay <p — > 
//(/ = u, d, c, s, t, b, e, /i, r, v e , v^, v T ), WW, ZZ, hh (for the radion decay width and it's branching 
ratio to several channels see Ref. [11]). 

Considering the NP contribution, the total decay width T(B S — ► p + p~) can be written as 



F(B S - M + /x-) = ^-|M TO t| 2 , (16) 



where |Mtot| 2 = \M$m + Mnp\ 2 (as given above (Eq. 14)). p c is the center of mass (c.o.m) 
momenta of the two charged muons in the B s rest frame and is given by 



(m| - ("V - "V-) 2 ) ( m | ~ ("V + "V) 2 ) 
^ c = 1 ' 2^ ' • (17) 

Finally the branching ratio BR(B S — > p + p~) is given by 

5i?(B s = ra s r(B a /x+/0, (18) 

where Tg (= 1.461 x 1CT 12 sec) is the life-time of the -B s meson [12]. 

3 Input parameters 

The decay amplitudes depend on the CKM matrix elements, wilson coefficients, quark and lepton 
masses and the non-perturbative input like decay constants. 

3.1 CKM matrix elements, quark masses, wilson coefficients and 
decay constant 

We adopt the Wolfenstein parametrization with parameters A, A, p and rj of the CKM matrix as 
below 



Vckm 



1 Vud Vus V ub 

V c d v cs V c b 
V V td V u V tb 



( 1-±A 2 A A\*(p-i v ) \ 

-A 1 - ±A 2 ^A 2 

V AX 3 (l-p-i7]) -AX 2 1 J 



(19) 
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We set A and A = sm6 c (9 C , the Cabibbo mixing angle) at 0.801 and 0.2265 respectively, in 
our analysis. Other relevant parameters are p = \/~p 2 + if cosj and 77 = \/p 2 + rf sinj, where 
\J~p 2 + if = 0.4048 and 7 ~ 62° [13]. The quark masses(in GeV unit) are being set at their current 
values i.e. m u = 0.2, = 0.2, m s = 0.2, m c = 1.4, mj, = 4.8, m< = 175 GeV [7] and = 0.105 
GeV. For the B s meson we use m-g = 5.369 GeV [12] and the decay constant fg = 0.20 GeV [14]. 
For our numerical calculation we use C10 = —4.5461 [7]. 

4 Numerical Analysis: Results and Discussions 

The branching ratio BR(B S — > p + p~)sM in the Standard Model is found to be two order smaller 
[6] than the present experimental upper bound 1.0 x 10~ 7 [5]. Assuming that the future data (not 
the upper bound) will still differ from the SM result, it is thus worthwhile to see whether radion 
can explain the discrepency. Our analysis is organized as follows: 

We first obtain the SM result. Using several numerical inputs (discussed in section 3), we find 

BR(B S ^ p + p-) SM = 3.17 x HT 9 , (20) 

which is two order smaller than the present experimental upper bound. Next we will see whether 
the radion can explain the above deficit or not. Taking into account the radion contribution along 
with the SM contribution, we find the total branching ratio which is given by Eq. (18) and is solely 
a function of the new physics parameters: (1) radion mass ms and (2) radion vev {(f)). In our 
analysis, we treat the radion as an ultra- light object whose mass ms varies in between 2m ^ GeV to 
mg, GeV and in this mass range it's decay width is found to be quite small. We set it at = 0.001 
GeV (see [11] for an elaborate discussion on radion decay width). 

Now let us define a quantity R = ^% Bs ^ l f M , - which is a function of ms and (6), the two 
NP parameters(as mentioned above). R, which is purely to be determined from the experiment, 
can be potentially quite significant in radion discovery. 

Since at present no data for BR(B S — ► fi + p~) is available (only upper bound exists), we take 
the conservative viewpoint and we use R to impose constraints on and (0). What we do is as 
follows: 

• We set R = 10, 50 and 100 and use those to obtain contour plots in the ms — (4>) plane, which 
are shown in Figure 4a. The topmost, middle and lowermost curves respectively stands for 
R = 10, 50 and 100. In Figure 4b we have plotted R as a function of m^ and {(j)). 

• From Figure 4a we can see that the lower bound on (0) for a given ms decreases with the 
increase in R. For R = 10(50, 100) and m & = 2 GeV, we find (0) > 918(564, 459) GeV. 
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Figure 4(a,b). Figure J^a corresponds to the contour plots in the mR{m$) -vR({cp)) plane. Here 
mR and vR are in GeV unit. The contour line starting from the top to bottom corresponds to 
R = 10, 50 and 100, respectively. In Figure 4b, we have shown the variation of R with mR and 
vR. 

• For a given R curve, we see that the lower bound on ((f)) increases with m^. For example, 
consider the R = 10 curve of Figure 4a. As m^ varies from 0.21 GeV to 5 GeV, ((f)) (the 
lower bound) changes from 844 GeV to 2620 GeV. Similarly for the R = 100 curve and for 
the above m^ range, the lower bound on ((f)) varies from 421 GeV to 1342 GeV. 

• In Figure 5 we have plotted BR(B S — > as a function of ((f)) for = 1 GeV (lower 
curve) and 4 GeV (upper curve), respectively. The horizontal line, the present experimental 
upper bound(= 1.0 x 10~ 7 ), suggest that the region allowed for ((f)) lies below the horizontal 
curve. This immediately allows one to obtain lower bound on (<p) and we find for = 1(4) 
GeV, the lower bound is about 605(945) GeV. 
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Figure 5. In Figure 5, the BR{B S — > ) as a function of vR(((f>)) (GeV) is shown. In y-axis 
we have actually plotted Log[BR(B s — > ^ + /i~)]. For example, BR = — 18(— 14) corresponds to 
1.523 x 10~ 8 (8.315 x 1(T 7 ). TTie lower(upper) curve corresponds to = 1(4) GeV, respectively. 
The horizontal curve (corresponding to BR(B S — > = 1.0 x 10~ 7 ,) represents the present 

experimental upper bound. The region below the horizontal curve is allowed for vR({<p) ). 

Invariant mass distribution 

The invariant mass distribution is an useful probe to see the NP signal and we now obtain 
such a distribution to see the radion signal. In Figure 6, we have plotted dT/dM^^ as a function of 
the M w (the di-muon invariant mass), respectively for ((f)) = 247, 500 GeV and 1 TeV (from top to 
bottom). The lowermost line corresponds to the SM background. Clearly for all {(/)) ranging from 
247 GeV to 1 TeV, the signal is way above the SM background and the height of the resonance 
peak increases. To see these consider the curves corresponding to = 1 GeV and ((f)) = 247, 500 
and 1000 GeV. One finds dT/dM^ as ~ 5.5 x 10" 13 , 8.55 x 10~ 13 and 1.03 x 10~ 12 , while for 
the SM one finds dT/dM^^ = 4.62 x 10~ 23 . The detectibility ofcourse depends on the clarity of 
the peak structure. Since the decay width goes as l/((p) 2 , with the increase of ((p), the resonance 
width decreases, one can have the better chance to see the resonance due to radion. We also 
see that with the increase of {(f)), resonance of higher radion mass gets lost and at ((f)) = 1 TeV, 
resonance corresponding to = 1, 2 GeV persists. We note that for (4>) = 2 TeV, the resonance 
corresponding to = 1 GeV only survives (which is not shown in the figure). 
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Figure 6. Invariant mass distribution dT/dM^ for the ultra-light radion signal with = 1,2,3 
and 4 GeV (from left to right for a given curve), where the lowermost curve corresponds to the SM 
background. In y-axis we have actually plotted Log[dT(B s — > fj, + fi~)/dM fJlfl ]. The topmost,middle 
and lowermost curves corresponds to (<fi) = 247, 500 GeV and 1 TeV, respectively. 
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Comparison with other studies 



We next compare our bound on (0) with that obtained from the muon anomaly data. The 
present BNL muon anomaly data tells a^ xpt ^ - a^ M ^ = (28 ± 10.5) x 1CT 10 . Fitting it with 
Sa^ (radion contribution to muon anomaly) the author in [15] found the following bound for a 
relatively light radion: for = 2 GeV, the lower bound on {4>) is ~ 500 GeV corresponding to 
da^ = 3.85 x 10~ 9 and ~ 760 GeV for <5a* = 1.75 x 10~ 9 . The upper bound on BR(B S -► n + n~) 
suggests a lower bound about 918 GeV (corresponding to R = 10) ~ 459 GeV (corresponding to 
R = 100) for ?7i0 = 2 GeV and is in agreement with that obtained from the muon anomaly data as 
mentioned above. 



5 Summary and Conclusion 

We analyse the B s — > fJ, + fJ,~ decay in the Randall-Sundrum model in the presence of an ultra-light 
stabilized radion. Using the present experimental bound of the branching ratio BR(B S — > 
(two order larger than the standard model result), we obtain contour plots in the — ((/)} plane. 
For BR(B S — > n + /J.~) = 1.0 x 10 -7 and radion mass = 0.21(5) GeV, we obtain the lower bound 
on ((f)) about 844(2620) GeV corresponding to R = 10 and about 421(1342) GeV corresponding 
to R = 100. We have shown the invariant mass (AL) distribution plot which clearly shows 
the existence of the radion signal. The SM background on which the radion signal lies is highly 
suppressed. Finally, we compare our result with the bound obtained from the BNL muon anomaly 
data and found some agreement between the two. 
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A SM amplitude of the process B s —> 

In this appendix, we calculate the square of the SM amplitude of B s — > fi + p~ (at the quark level 
it yields bs — ► fi + p~ process) of Eq.(9). We define p b , p s , pi and p 2 to be the momenta of the 
b-quark, s-quark, fi + and pT , respectively and m^, the muon mass and q = p b + Ps = Pi +P2- The 
individual amplitude-square elements are listed below: 



l-Miol 2 = C'\C w \ 2 [Am^m b 2 - 4{ Pl .p 2 )m b 2 + 4m/m s 2 - 4{pi.p 2 )m s 2 

+8(Pi-Pb)(P2-Pb) + 8{p 1 .p s )(p 2 .pb) + 8(p 1 .p b ){p 2 .p s ) + 8(pi.p s )(p 2 .p s ) 
+8m^ 2 (p b . Ps ) - 8( Pl .p 2 )(p b .p s )}. (21) 

where the prefactor C is defined as follows 

C' = (^^) 2 \V tb \ 2 \V ts \ 2 . (22) 
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B NP contribution 



B.l Interaction vertex of Radion and SM particle 




G ^ 

'' >k 2 2i_[k k - 2^M 2 1 

^ k i 12 b « J 

G" " 




Y [ 1 - Y 5 ] Vj, 



w, 



2i 



x v>k <(b> 

v G Y 



M w [ k G - k w ] 




1 9 [mj ( l- Y 5 ) -mj ( i + Y 5 )] Vj , 

2/2 M w 



Figure 1: The relevant diagrams of the interaction vertex of Radion and SM particles in the 
gauge. 
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B.2 Feynman diagrams of the Radion penguin with SM particles. 




b r 

-> C 



■A 



->- 



1 



-> f- 



\ S 

— ' >- 



(g) 



(h) 



(i) 



b / 



(j) 



Figure 2: The Feynman diagrams of the Radion penguin with SM particles in the general 
Rt: gauge . 
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B.3 NP amplitude of the process B s — > fi + fi 

In this appendix, we calculate the square of the NP amplitude of B s — > ^ + [i>~ of Eq.(14). The 
individual amplitude-squared elements are given by 



where 



In above 



Cl(j>i-P2 ~ rn^ 2 ), 



2Re{M\ M N p) = 0, 



= \V tb \'\V ts 
m 2 rl 



1 



^ v^(0) J \m h + m s ) ( (? 2_ m 2 )2 + r 2 m 2 ) 



27T 2 (</))yo 



F{x)dx. 



F{x) = 4(Log(A 2 /A 2 )-l) + ^ + 

A = 4:ir((f>}, A 2 = xm 2 + (1 - x)M w 2 . 

We worked in the 't-Hooft Feynman gauge (£ = 1). 



mi( 2 (1 — a;) Mi 



(23) 
(24) 
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